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ABSTRACT 


A telephone  channel  simulator  has  been  implemented  on  the  Lincoln 
Digital  Voice  Terminal  (LDVT) . This  technical  note  first  reviews  the 
various  distortions  that  occur  in  telephone  channels  and  then  describes 
the  implementation.  This  work  was  undertaken  at  the  suggestion  of 
Mr.  Ronald  Sonderegger  of  the  Defense  Communications  Engineering  Center 
who  observed  that  a real-time  simulation  on  the  LDVT  should  be  far 
less  costly  than  previous  implementations.  The  results  obtained  indicate 
that  the  potential  of  the  LDVT  goes  beyond  the  original  intent  of  its 
use  as  a speech  digitizer. 
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INTRODUCTION 


There  has  been  growing  interest  recently  in  the  feasibility  of  inter- 
connecting the  telephone  system  with  the  various  digital  vocoders  that  are 
currently  being  developed.  The  obvious  possible  uses  of  the  telephone  system 
are:  1)  to  transmit  speech  from  a local  source  to  the  input  of  a vocoder,  and 
2)  to  transmit  the  bits  from  a vocoder  analyzer  to  a vocoder  synthesizer. 

Each  of  these  possible  functions  has  a set  of  problems  associated  with  it,  due 
to  the  various  distortions  and  interferences  in  a telephone  line.  For  example, 
it  is  generally  known  that  pitch  detection  algorithms,  upon  which  many  vocoders 
depend,  break  down  when  the  input  is  telephone  speech.  When  bits  are 
transmitted  over  a telephone  line,  they  must  be  converted  to  a type 
of  signal  which  can  be  sent  with  a minimal  number  of  errors.  This  task  is 
assumed  by  a modem,  of  which  there  are  many  varieties  on  the  market.  What 
makes  this  task  difficult  is  the  fact  that  the  telephone  channel  makes  no 
attempt  to  correct  for  phase  distortion,  since  the  human  ear  is  presumed  to 
be  insensitive  to  phase  delay. 

By  simulating  the  telephone  channel's  characteristic  distortion  and  inter- 
ference in  real  time  in  a digital  computer,  it  should  be  possible  to  investigate 
the  reasons  for  the  problems  encountered  when  actual  telephone  lines  are  used, 
both  as  inputs  to  a vocoder  and  as  transmission  media  for  bits.  Using  the 
simulation,  one  can  hopefully  develop  more  robust  modems  and  vocoders,  and 
one  can  determine  how  much  of  a given  type  of  interference  is  required  before 
a given  system  breaks  down. 
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This  report  begins  with  a discussion  of  how  the  telephone  system  is  con- 
structed. Each  of  the  various  types  of  interference,  Gaussian  noise,  frequency 
attenuation  and  phase  delay,  quadrature  distortion,  nonlinear  distortion,  echo, 
and  cross  talk,  is  then  discussed  in  some  detail.  Following  that,  a description 
is  given  of  how  the  measurements  were  obtained  for  the  simulation.  Finally, 
the  details  of  the  actual  simulation  are  described  fully,  including  an  indication 
of  those  areas  in  which  the  simulation  is  a simplification  of  the  actual  situation. 

I . Overview-Description  of  a Telephone  Channel 

In  a typical  telephone  channel1,  the  signal  is  first  filtered,  then  mod- 
ulated up  to  some  carrier  frequency,  transmitted  with  multiplexing  through 
a series  of  cables  and  repeaters  and  finally  filtered  and  demodulated  at  the 
receiver.  Various  distortions  are  introduced  into  the  signal  at  all  parts  of 
the  system.  The  filters  pass  only  energy  between  300  and  3000  Hz,  and  introduce 
both  amplitude  and  phase  distortion  in  the  pass  band.  The  modulation  demodulation 
process  introduces  quadrature  distortion  whenever  there  is  an  offset  between 
the  modulating  frequency  and  the  demodulating  frequency.  Cross  talk  is  a result 
of  both  frequency  and  space  division  multiplexing,  whenever  the  filtering  and/or 
shielding  is  inadequate.  In  travelling  down  the  cable,  the  signal  becomes 
attenuated  and  therefore  the  repeaters  are  necessary  to  amplify  it  at  various 
points  in  the  transmission  lines.  The  repeaters  are  a source  of  both  Gaussian 
noise  and  nonlinear  distortion.  When  a signal  is  switched  from  one  carrier 
to  another  there  is  introduced  a sudden  change  in  the  phase  relationship,  as 
well  as  transient  (impulse)  noise.  Another  source  of  impulse  noise  is  lightning 
and  corona  type  discharges. 
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II.  SSB  Modulation 


Typically,  voice  signals  are  frequency  multiplexed  over  a telephone 
line  using  single  sideband  suppressed  carrier  (SSBSC)  amplitude  modulation,2’3 
as  shown  in  Figure  1.  In  order  to  set  the  stage  for  our  subsequent  discussion 
of  phase  distortion,  SSBSC  is  reviewed  in  this  section. 

The  signal  is  first  bandpass  filtered  from  300  to  3000  Hz  and  then  modulated 
up  by  a carrier  cosine  whose  frequency  is  a multiple  of  4000  Hz.  The 
resulting  signal  is  then  bandpass  filtered  to  eliminate  any  energy  below 
Since  the  original  signal  was  real,  this  entails  no  loss  of  information.  The 
signal  is  then  frequency  multiplexed  with  several  other  signals  at  different 
carrier  frequencies  (also  multiples  of  4 kHz)  and  is  detected  at  the  receiver 
by  means  of  a bandpass  filter  from  +_  u)c  to  ^(u>c+4000  Hz),  followed  by  demodulation 
using  a cosine  also  at  frequency  However,  no  attempt  is  made  to  assure  that  the 

modulating  and  demodulating  cosine  waves  are  in  phase  with  each  other . In  fact  the  mod- 
ulating and  demodulating  frequencies  are  rarely  exactly  the  sam^,  and  any  slight  diff- 
erence in  frequency  will  cause  a gradual  drift  into  and  out  of  phase. 

When  the  modulating  and  demodulating  cosines  are  exactly  90°  out  of  phase 
the  resulting  signal  which  is  received  is  the  quadrature  component,  or  Hilbert 
Transform,1*  of  the  original  signal.  The  quadrature  component  can  be  understood 
from  several  different  viewpoints.  If  one  represents  the  original  signal  f(t) 
in  terms  of  its  Fourier  expansion,  i.e.,  as  a sum  of  cosines  at  different 
frequencies  with  differing  amnlitudes  and  phases,  then  the  Hilbert  Transform 

/v 

f(t)  could  be  constructed  by  converting  all  of  the  cosines  to  sines  and  adding 
them  up.  An  impulse,  for  example,  is  a sum  of  cosines  all  in  phase  at  9 = 0, 
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at  that  point  in  time  when  the  impulse  occurred.  If  all  of  these  cosines  were 
converted  to  sines,  every  function  would  be  zero  at  the  time  of  the  impulse, 
all  would  be  negative  immediately  before  it,  and  all  would  be  positive 
immediately  after.  Some  distance  away  on  either  side,  the  sines  would  be  at 
random  phase  with  respect  to  one  another,  and  hence  the  net  sum  at  any 
point  would  be  0.  The  result  is  the  function  shown  in  Figure  2. 

One  way  to  acquire  the  Hilbert  Transform  is  to  convolve  the  signal  with 
the  function  of  Figure  2.  Another  way  is  to  take  the  Fourier  transform, 
multiply  positive  frequency  by  -j  and  negative  frequency  by  +j  and  take  the 
inverse  Fourier  transform.  It  is  clear  that  such  a filter  would  convert 
a cosine  at  any  frequency,  to,  into  a sine,  since  the  transform  of  a cosine 
is  two  positive  real  impulses  at  j^to  and  the  transform  of  a sine  is  a negative 
imaginary  impulse  at  + to , and  a positive  imaginary  impulse  at  -urffe^-e  •*W)/2j  = 
sin  to) . 

Supposing  one  had  available  a complex  signal  z(t)  whose  spectrum  was 
zero  for  negative  frequency  and  identical  to  the  spectrum  of  a given  real  signal 
f(t)  for  positive  frequency.  According  to  the  previous  discussion,  such  a 
function  could  be  generated  as  follows: 

z(t)  = 1/2  (f(t)  + j f(t))  (Figure  3) 

If  one  were  to  multiply  z(t)  by  e^Wct,  the  transform  of  the  resulting 
function  would  be  simply  Z(<o)  translated  up  by  a frequency  to^,  i.e.,  Z(w-w  ). 
Now,  were  one  to  take  the  real  part  of  this  signal  and  determine  its  transform 
the  result  would  be  the  even  part  of  Zfoo-o^)  which  (Figure  4)  is  clearly  the 
same  as  what  would  be  obtained  by  multiplying  f(t)  by  cos  and  filtering  out 
energy  below  oo^.  Remembering  that  the  even  part  of  a transform  is  the  transform 
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of  the  real  part  of  the  function,  we  have: 
fs(t)  = 2 Re  [z(t)  ejatt] 

= Re  [(f(t)  + j f(t))  e^) 

= [f(t)  cos  a)ct  - f(t)  sin  u>  t] . 

Now  it  should  be  straightforward  to  explain  the  effect  of  a phase  difference 
between  the  modulating  and  demodulating  carriers. 

fQ(t)  = [f(t)  cos  to ct  - f(t)  sinu^t)]  cos  (u^t  + <j> ) (Low  Pass) 

A 

= f(t)  cos  wct  cos  (u^t  + <f> ) - f(t)  sin  o^t  cos  (o^t  + <}> ) (Low  Pass) 

= f(t)  cos  w^t  cos  oi^t  cos  <t>  - f(t)  cos  to^t  sin  w^t  sin  <p 

A A 

- f(t)  sin  w^t  cos  oi^t  cos  (|)  + f(t)  sin  co^t  sin  to ^t  sin  <p  (Low  Pass) 

Recalling  that 

cos2  6 = 1/2  (1  + cos  2 6) 

sin2  0 = 1/2  (1  - cos  2 0) 

sin  0 cos  0=  1/2  (sin  2 0) 

and  that  the  sin  2 0,  cos  2 0 terms  will  all  be  clobbered  by  low  pass  filtering 
at  the  output,  we  have  finally 

fQ(t)  = 1/2  (f(t)  cos  <}>+  f(t)  sin  4> ) 

When  <(>=  90°,  fQ(t)  = 1/2  f(t),  which  is  to  say  that  only  the  quadrature 
term  is  received. 

Ill . Disturbances  in  the  Transmission  System 
A.  Quadrature  Distortion 

Armed  with  a knowledge  of  the  mathematics  of  SSB  carrier  systems,  it 
is  now  simple  to  understand  the  various  phase  distortions  that  occur  in  the 
telephone  systems.  If  there  is  a frequency  difference , between  the  modulating 
and  demodulating  waves,  we  would  have: 
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fQ(t)  = (f(t)  cos  wct  - f(t)  sin  W£t)  cos  ((wc  + Aw)t)  (Low  Pass) 

The  contribution  to  the  phase  difference,  <j)  , by  the  frequency  offset, 
u , is  a function  of  time: 

<J>  Q = ( Aw  ) t 

which  accounts  for  a slow  drift  into  and  out  of  phase. 

Phase  jitter  is  the  term  used  to  describe  the  interference  of  a sinusoidal 
noise  (frequently  60  cycles)  in  the  phase  of  a carrier.  Sixty  cycle  inter- 
ference shows  up  everywhere  in  the  system  but  in  general  is  harmless  because 
of  the  300  Hz  lower  limit  of  the  filters.  However,  60  cycle  interference 
in  the  generation  of  a cosine  causes  a phase  distortion  of  the  cosine  waveform 
which  appears  as  a very  low  frequency  modulation  (Figure  5) . The  phase 
jitter  introduces  an  additive  component  to  the  phase  which  can  be  expressed 
as 


d>  . = A.  cos  a)  .t 
3 3 3 

where  is  the  amplitude  of  the  jitter  in  degrees  anduij  is  usually  60  cycles 
(50  cycles  in  Europe) . 

The  final  phase  effect  is  phase  hits,  or  sudden  large  shifts  in  phase, 
which  occur  as  a result  of  switching  two  carrier  supplies  not  in  phase.  Some- 
times a phase  coherent  detector  corrects  the  situation  some  time  later.  On 
other  occasions  the  original  phase  relationship  is  never  restored. 

The  net  phase  relationship  is  a combination  of  the  three  effects: 

<f>  = <f> Q + <f>  j + <f>h 
and  the  final  output  signal  is 

A 

fQ(t)  = f(t)  cos  <p  + f(t)  sin  <J)  . 


6 


Filtering 


A major  part  of  the  distortion  in  a telephone  signal  has  to  do  with 
the  frequency  response  and  phase  delay  characteristics  of  the  filter  used 
to  assure  that  the  signal  remains  in  band.  Theoretically,  in  multiplexing 
at  4000  Hz  intervals,  one  could  pass  all  frequencies  from  0 to  4000  Hz 
and  avoid  cross  talk.  In  practice,  filters  never  cut  off  sharply,  and 
so  a healthy  compromise  is  to  pass  the  signal  from  about  300  to  3000  Hz. 

In  addition,  there  tends  to  be  considerable  noise  interference  at  low  frequencies 
(for  example,  60  cycle  hum)  so  that  it  is  advantageous  from  a signal  to  noise 
ratio  argument  to  remove  low  frequencies.  The  characteristics  of  the  filter 
in  the  pass  band  represent  a tradeoff  between  cost  and  response.  In  particular 
not  much  attention  was  paid  to  phase  delay  distortion,  as  the  ear  is  insensitive 
to  phase. 

C.  Nonlinear  Distortion 

In  the  process  of  travelling  down  a coaxial  cable,  a signal's  energy 
is  gradually  attenuated  due  to  losses  in  the  transmission  line  and 
therefore  it  is  necessary  to  amplify  the  signal  at  certain  points  in  the 
transmission  path  to  assure  an  adequate  level  at  the  receiver.  The  devices 
which  achieve  the  amplification  are  called  repeaters,  and  are  generally  made 
up  of  two-port  amplifiers.  The  voltage  transfer  characteristics  of  a 
generalized  two-port  are  shown  in  Figure  6.  According  to  the  Fourier 
series  expansion,  the  transfer  characteristics  can  be  expressed  as: 


where  a^  is  the  gain  and  the  other  a's  are  the  nonlinear  distortion  coefficients. 


3 

+ a_e.  + 
3 i 
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The  sources  of  the  nonlinear  distortions  are  several,  and  much  has  been 
written  on  the  subject  of  repeater  design  (see,  for  example.  Bell,  pp.  396-421, 
Ref.  1) . Transistors  have  an  inherent  nonlinearity  dependent  on  the  source 
and  load  impedances  and  the  selected  d-c  bias  conditions.  Negative  feedback 
in  the  amplifier  complicates  the  nonlinear  relationship.  The  nonlinearities 
in  a repeater  are,  in  general,  frequency  dependent  and  amplitude  dependent, 
and  the  calculation  of  the  precise  response  for  a given  system  design  is  quite 
complicated. 

D.  Gaussian  and  Impulse  Noise 

Gaussian  noise  is  an  unavoidable  side  product  of  networks  and  devices. 

Two  common  types  of  noise  in  a circuit  are  thermal  noise  and  shot  noise.5 
Arguing  from  the  central  limit  theorem,  both  can  be  assumed  to  be  Gaussian  and 
white  at  the  source.  Thermal  noise  is  present  in  any  conductor,  for  example, 
a resistor,  and  is  due  to  the  thermal  interaction  between  free  electrons  and 
vibrating  ions.  Its  available  power  is  directly  proportional  to  the  product 

of  the  bandwidth  of  the  system  and  the  temperature  of  the  source.  Shot  noise 

v 

is  due  to  the  discrete  nature  of  electron  flow  and  is  present  in  most  active 
devices  (transistors  and  diodes) . Its  amplitude  is  proportional  to  the  square 
root  of  the  current  and  therefore  is  dependent  on  signal  level. 

The  components  of  the  repeaters  are  a major  source  of  the  Gaussian  noise. 
By  the  time  the  noise  reaches  the  receiver,  it  is  no  longer  white,  because 
it  has  been  shaped  by  the  filters  at  demodulation. 

In  addition  to  Gaussian  noise  there  is  occasionally  a burst  of  noise 
on  a line  whose  amplitude  far  exceeds  the  average  noise  level  of  the  system. 
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Such  noise  is  referred  to  as  "impulse  hits"  and  is  generally  caused  by 
switching  transients  in  central  offices  or  from  corona  type  discharges 
(electrical  discharges  in  the  air  surrounding  a high  potential  line)  that 
occur  along  a repeated  line. 

E.  Echo  and  Cross  Talk  * 

To  be  complete  in  a report  on  telephone  channel  characteristics  one 
would  have  to  include  at  least  a brief  discussion  of  echo  and  cross  talk,  even 
though  these  two  effects  were  not  included  in  the  digital  simulation.  An 
echo  may  be  produced  whenever  there  is  an  impedance  discontinuity.  It  is 
most  commonly  caused  by  a return  of  a talker's  signal  through  the  channel 
to  which  he  is  listening.  For  short  distances  the  effect  is  indistinguishable 
from  side  tone  and  is  therefore  not  annoying.  Echo  suppressor  circuitry  has 
been  added  to  long  distance  lines  to  attenuate  the  signal  on  the  return 
path  when  the  received  signal  amplitude  is  high.  The  result  is  that  when 
both  people  talk  only  one  is  heard,  and  sometimes  the  beginnings  of  words  are 
clipped.  There  is,  however,  a "tone  disabler"  section  of  the  echo  suppressor 
circuit  which  permits  the  mechanism  to  be  shut  off  when  data  is  being  transmitted. 

Cross  talk  is  caused  by  coupling  losses  between  two  active  circuits. 
Transmittance  cross  talk  occurs  because  of  inadequate  design  of  modulators 
and  filters  in  frequency  multiplexing.  Coupling  cross  talk  is  caused  by 
electromagnetic  coupling  between  two  physically  isolated  circuits.  The 
result,  of  course,  is  that  the  listener  hears  another  conversation  in  the 
background,  which  may  or  may  not  be  intelligible. 
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IV.  Measuring  the  Telephone  Channel 


Before  the  telephone  channel  could  be  simulated,  it  was  necessary  to  make 
extensive  measurements  on  thousands  of  lines  in  order  to  determine  concrete 
values  to  be  used  in  the  various  components  of  the  simulation.  From  the 
study  data  set,  histograms  were  compiled  for  all  of  the  various  parameters, 
such  as  frequency  offset  and  level  of  Gaussian  noise,  for  Continental  U.S. 
(Conus)  and  European  voice  and  data  grade  lines.  Table  1 is  a chart  of  the 
resulting  numbers  for  the  various  degradations  (excepting  filter  frequency 
response)  .6  The  term  "mid"  is  used  to  refer  to  the  50th  percentile  on  the 
histogram  and  "poor"  to  refer  to  the  90th  percentile;  i.e.,  90%  of  the  Conus 
voice  grade  channels  measured  were  better  than  "Conus  poor  voice". 

Gaussian  noise  and  nonlinear  distortion  were  measured  in  units  of  dBmc, 

meaning  dB's  relative  to  1 milliwatt,  using  a special  C-message  weighting 

curve ’(Figure  7).  This  frequency  weighting  curve  was  determined  empirically 

by  means  of  several  subjective  listener  tests  and  is  intended  to  reflect 

the  amount  of  listener  annoyance  for  noises  at  different  frequencies.  The 

unit  of  measurement  used  for  impulse  noise  is  dBrnc,  where  rn  stands  for 

-12 

reference  noise,  which  is  -90  dBm,  or  10  watts  at  1000  Hz. 

The  device  used  by  the  telephone  company  to  measure  Gaussian  noise  has 
built  in,  in  addition  to  the  C-message  weighting  curve,  a mechanism  for 
attenuating  impulse  noise  of  a sufficiently  short  duration.  The  argument  again 
is  one  of  ear  sensitivity,  as  the  human  ear  does  not  fully  perceive  the 
power  in  noises  of  duration  less  than  200  msec. 

The  telephone  company  has  available  another  device,  called  an  impulse 
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counter,  for  the  purposes  of  measuring  the  amount  of  impulse  noise  present,  as 
such  noise  is  far  more  destructive  than  Gaussian  noise  in  introducing  errors 
in  the  transmission  of  bits  through  a modem  and  a channel.  The  counter  consists 
of  a weighting  network,  a rectifier,  a threshold  detector,  and  a counter  of 
events  above  threshold.  Several  impulse  counters  can  be  used  simultaneously  at 
different  thresholds  to  obtain  information  about  the  distribution  of  the 
magnitudes  of  the  impulses. 

To  measure  the  nonlinear  distortion,9  a pure  tone  at  a fixed  level  is 
transmitted,  and  the  amount  of  energy  received  at  the  2nd  and  3rd  harmonics 
is  measured.  From  these  values,  one  can  deduce  the  coefficients  a,  and  a^ 
of  the  squared  and  cubed  terms : 

cos^O  = 1/2  (1  + cos  2 9) 

3 

cos  0 = cos0  cos  20  = 1/2  (cos0  + cos0  cos  20  ) 

= 1/2  cos0  + 1/4  (cos  30  + cos0  ) 

Because  of  this  relationship  between  nonlinearities  and  harmonics,  nonlinear 
distortion  is  often  referred  to  as  harmonic  distortion. 

V.  Simplifications  Used  in  the  Simulation 

The  parameters  of  the  simulation  were  set  so  as  to  emulate,  as  closely 
as  possible,  the  results  found  from  the  statistical  study.  In  the  case  of 
certain  measures,  such  as  frequency  offset  and  phase  jitter,  this  was  a straight- 
forward process.  However,  in  the  simulation  of  Gaussian  noise,  impulse  hits, 
phase  hits,  and  nonlinear  distortion,  certain  (sometimes  gross)  simplifications 
were  used. 
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The  numbers  used  for  Gaussian  noise,  impulse  noise,  and  nonlinear 
distortion  in  the  simulation  were  determined  empirically  by  adjusting  constants 
until  the  desired  value  was  obtained  at  the  output  by  the  appropriate  measuring 
device  (described  in  section  IV). 

Gaussian  white  noise  generated  at  various  devices  in  the  transmission  line 
is  no  longer  white  by  the  time  it  reaches  the  receiver  because  of  the  filter 
transfer  characteristics.  In  the  simulation,  Gaussian  noise  is  added  at  the 
input,  so  that  it  will  be  shaped  by  the  filter  of  the  system.  The  level  of  the 
Gaussian  noise  in  the  telephone  system  is  dependent  upon  the  amplitude  of 
the  input  signal,  whereas  in  the  simulation  noise  is  kept  at  a fixed  level. 

Impulse  hits  and  phase  hits  are  simulated  to  occur  at  fixed  intervals 
and  at  fixed  duration  and  level,  so  as  to  correspond  in  frequency  of  occurrence 
and  average  amplitude  to  the  given  telephone  channel  characteristics. 

Nonlinear  distortion  is  simulated  by  squaring  and  cubing  the  final  output 
of  the  system.  The  coefficients  a^  and  a^  remain  fixed  and  are  equal,  even 
though  in  the  actual  telephone  system  these  two  coefficients  would  in  general 
be  dependent  on  both  amplitude  and  frequency.  The  appropriate  gain  for  the 
squared  and  cubed  terms  was  of  course  determined  so  as  to  match  the  desired 
meter  measurement.  This  was  done  by  sending  a pure  tone  at  700  Hz  through  the 
digital  distorter  and  measuring  the  output  of  a notched  filter  (to  remove 
700  Hz)  in  dBrnc. 

VI . Implementation  on  the  LDVT 

The  Lincoln  Digital  Voice  Terminal  (LDVT)  is  a small  high  speed  computer 
which  has  proven  capable  of  realizing  in  real-time  a variety  of  algorithms 


12 


for  low  bit  rate  digital  speech  transmission  1 °( from  2400  bits/sec  to  16000 
bits/sec).  The  computer  consists  of  a 1024  word  program  memory,  Mp,  with  a 
cycle  time  of  55  nsec,  a 512  word  data  memory,  Md,  and  a 2048  word  outboard 
memory,  Mx,  an  input/output  device  from  which  data  is  accessible  in  a few 
instruction  cycles.  The  LDVT  has  a very  simple  instruction  set  and  is  quite 
convenient  to  program. 

The  telephone  channel  simulator  represents  the  first  attempt  to  use  the 
LDVT  for  something  other  than  a vocoder.  A block  diagram  of  the  complete 
simualtion  is  given  in  Figure  8.  The  input  speech  is  low  pass  filtered  to 
remove  energy  above  5 kHz  and  sampled  at  10  kHz  using  a 12  bit  A/D  converter. 
Impulse  and  Gaussian  noise  are  generated  digitally  and  added  to  the  input  samples 

The  frequency  response  and  delay  characteristics  of  the  telephone  filter 
and  the  Hilbert  transform  to  obtain  the  quadrature  component  are  both  realized 
by  means  of  high  speed  convolution  techniques.11  After  a 256  point  FFT  of  the 
input  samples  (256  points  is  adequate  as  long  as  the  combined  impulse  response 
of  the  filter  and  Hilbert  transform  is  less  than  12.8  msec),  negative  frequency 
is  zeroed  out  and  positive  frequency  is  complex  multiplied  by  the  FFT  of  the 
desired  impulse  response.  A 256  point  IFFT  now  yields  the  filtered  speech  f(t) 

A 

in  the  real  data  buffer  and  the  filtered  quadrature  component,  f(t),  in  the 
imaginary  data  buffer  (refer  to  the  previous  discussion  on  SSB  modulation) . 

Only  the  second  half  of  the  data  is  good,  because  of  circular  convolution, 
and  therefore  an  overlap  of  128  samples  is  necessary  between  frames. 

The  phase, $ , is  computed  as  the  sum  of  the  contributions  of  phase  jitter, 
frequency  offset,  and  phase  hits.  The  output  of  the  "modulator  demodulator" 
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complex  is  then  computed  as: 

g(n)  = f(n)  cos  tj>  + f(n)  sin  ij>  . 

The  final  step  in  the  simulation  is  to  subtract  from  g(n)  the  nonlinear 
distortion  term: 

s~(n)  = g (n)  - D(g^(n)  + g3(n)). 

There  is  a special  buffer  in  Md  (data  memory)  which  is  set  aside  as  the 
parameters  of  the  system  (excepting  filter  frequency  response)  and  which  can 
be  modified  under  user  control.  Figure  9 is  a copy  of  that  section  of  the 
program,  including  the  appropriate  settings  for  Conus  poor  voice.  The  same 
mnemonics  are  used  in  Figure  8,  where  it  should  be  clear  how  the  various  para- 
meters are  being  used. 

Table  2 gives  the  values  used  for  the  various  parameters  in  the  simula- 
tion of  the  eight  channels.  The  values  for  the  parameters  FRQOFF , PKTOPK, 
JITTER,  PHSHIT,  PKTFRC  and  OFFFRC  are  simply  copied  over  from  Table  1.  The 
values  for  PHSRTE  and  IMPRTE  are  determined  as  the  reciprocal  of  #/15  min, 
and  converted  from  units  of  minutes  to  units  of  .01  sec.  The  parameters  PHSDUR 
and  IMPDUR  are  left  unspecified  in  Table  1,  and  are  set  arbitrarily  to  .01 
sec  in  Table  2 for  all  channels.  The  parameters  GAUSML,  HARM,  and  IMPAMP  are 
dimensionless  fractions  determined  empirically  so  as  to  read  the  correct  meter 
reading  at  the  output  in  units  of  dBmc  or  dBrnc. 

The  only  other  information  that  is  needed  to  simulate  each  of  the  channels 
is  the  frequency  domain  characteristics  of  the  telephone  filter.  Tables 
were  given  for  each  channel  of  256  real  and  256  imaginary  frequency  components 
of  the  desired  filter,  spaced  by  20  Hz,  spanning  0 to  5 kHz.  Since  the  DVT 
implementation  used  a 256  rather  than  512  point  FFT,  alternate  samples 
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were  ignored  in  entering  the  tables  into  the  computer.  It  should  be  noted  here 
that  it  would  have  been  essentially  impossible  to  implement  a 512  point  FFT 
in  the  DVT  due  to  its  current  memory  limitations,  but  that  the  12.8  msec  window 
appears  to  be  adequate  for  the  filters  simulated. 

The  main  body  of  the  program  is  the  FFT-IFFT  computation  which 
requires  several  buffers  in  Mx  (outboard  memory),  two  128  point  buffers  in 
Md  and  a complex  interchange  of  data  between  Mx  and  Md,  The  program 

takes  advantage  of  certain  characteristics  of  the  data  to  reduce  both  time  and 
memory  requirements.  By  keeping  the  FFT  size  down  to  128,  one  avoids  the  necessity 
of  referencing  Mx  in  the  inner  loop  which  would  greatly  increase  the  time 
required.  Therefore  the  forward  FFT  is  realized  (since  the  input  data 
is  real)  by  packing  even  numbered  samples  in  the  real  Md  buffer, 

MDREAL,  and  odd  numbered  samples  in  the  imaginary  Md  buffer,  MDIMAG, 
and  doing  a 128  point  FFT  followed  by  even  odd  separation  and  a final  stage 
of  the  256  point  FFT.  The  first  stage  of  the  inverse  FFT  is  nonexistent,  since 
the  second  half  of  the  data  is  zero,  and  therefore  it  can  be  conveniently 
split  into  two  128  point  FFT's,  with  the  first  one  yielding  the  even  numbered 
output  samples  and  the  second  one  (by  beginning  each  stage  with  the  coefficient 
index  set  at  half  its  increment  instead  of  at  zero)  yielding  the  odd  numbered 
output  samples. 

Because  of  this  even  numbered  odd  numbered  sorting  at  both  the  output 
and  the  input,  it  is  convenient  to  store  all  buffers  of  speech  in  Mx  in  this 
peculiar  fashion  of  even  numbered  samples  in  one  buffer  and  odd  numbered  samples 
in  another.  The  buffers  required  in  Mx  are  listed  in  Table  2.  The  even  num- 
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bered  input  samples  are  stored  in  EVIO  and  the  odd  numbered  in  ODIO.  At 
each  A/D  interrupt  the  next  sample  is  fetched  out  of  either  EVIO  or  ODIO  alter- 
nately and  sent  to  the  D/A  converters  and  the  new  sample  from  the  A/D  converter 
is  written  over  the  same  place  in  Mx.  At  the  beginning  of  a new  frame,  the 
most  recent  128  samples  from  the  A/D  converter  have  filled  up  EVIO  and  ODIO, 
and  so  these  buffers  now  become  EVNEW  and  ODNEK  respectively.  At  the  same 
time,  the  former  EVNEW,  ODNEW  become  EVOLD,  ODOLD,  and  the  former  EVOLD,  ODOLD 

A 

whose  contents  had  been  written  over  by  s(n)  in  the  course  of  the  previous 
frame,  become  EVIO  and  ODIO,  ready  to  be  sent  to  the  D/A  converter. 

The  first  step  in  a new  frame  is  to  fill  MDREAL  with  128  even  numbered 
new  samples  from  EVOLD/EVNEW  and  to  fill  MDIMAG  with  the  odd  numbered 
samples  from  ODOLD/ODNEW.  Now  a 128  point  normal  to  bit  reversed  decimation 
in  frequency  FFT  followed  by  bit  reversed  even  odd  separation  and 
a final  stage  of  the  256  point  FFT  can  all  be  realized  in  place  in  Md.  The 
resulting  data  (the  first  128  samples  of  the  256  point  FFT)  is  then  complex 
multiplied  by  the  bit  reversed  filter  coefficients  (which  have  been  stored  in 
Mx) . The  resulting  filtered  spectrum  is  then  saved  in  Mx:  the  real  data  is 
stored  over  ODOLD/EVOLD  and  the  imaginary  data  is  saved  in  a 128  point  buffer, 
MX I MAG.  A 128  point  IFFT  at  this  point,  of  the  data  in  Md,  yields  the 

A 

even  numbered  f(n)  samples  in  MDREAL  and  the  even  numbered  f(n)  (quadrature 
component)  samples  in  MDIMAG. 

Only  the  2nd  half  of  the  data  is  good,  because  of  circular  convolution, 

A 

and  these  64  even  numbered  f(n)  and  even  numbered  f(n)  are  saved  in  Mx.  Mean- 
while, the  original  real  and  imaginary  outputs  of  the  forward  FFT  filtered 
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spectrum  are  fetched  back  from  Mx  and  another  128  point  IFFT,  beginning  each 
stage  with  the  coefficient  index  set  at  half  its  increment  instead  of  at  0, 

A 

yields  the  odd  numbered  samples  of  f(n)  and  of  f(n)  in  Md.  Again,  the  first  64 

A 

of  each  are  garbage.  The  64  even  numbered  samples  of  f(n)  and  of  f(n)  can  now 
be  fetched  back  from  Mx  and  stored  over  the  64  garbage  points  in  each  of 

MDREAL  and  MDIMAG.  Now  all  of  the  pertinent  data  resides  in  Md  and  all 

that  is  left  is  the  introduction  of  quadrature  distortion  and  nonlinear 
distortion: 

g(n)  = f(n)  cos  <J>  + f(n)  sin  <J> 

s (n)  = g(n)  - D (g2  (n)  + g3  (n)) 

A 

This  can  conveniently  be  done  at  this  time  and  the  resulting  s (n)  evens  can 

A 

be  stored  over  EVOLD  and  s(n),  odds,  over  ODOLD. 

When  the  interrupt  routine  has  been  serviced  128  times,  resulting  in  64 
new  samples  in  each  of  EVIO,  and  ODIO,  it  is  time  to  swap  buffer  pointers  again, 

leaving  s(n)  evens  in  EVIO  and  s(n)  odds  in  ODIO,  as  desired, and  we  have 
come  full  circle. 

The  assurance  of  adequate  accuracy  without  overflow  in  the  FFT-IFFT 
complex  at  first  posed  some  difficulty,  as  a check  for  overflow  and  correction 
in  the  inner  loop  proves  very  costly  in  terms  of  time.  The  scheme  that  was 
finally  decided  upon  was  the  following: 

1.  Shift  the  input  buffer  of  the  forward  FFT  in  Md  as  far  left 
as  possible  as  a block  so  that  the  largest  sample  does  not  overflow. 

2.  Scale  down  the  data  as  a block  by  1/2  at  each  stage  of  the 
inner  loop  of  the  forward  FFT. 
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3.  Scale  the  output  data  in  Md  from  the  forward  FFT  as  far  left 

as  possible. 

N 

4a.  If  the  data  has  been  scaled  up  by  more  than  2 , scale  it  down 

by  1/2  at  each  stage  of  the  IFFT  and  scale  it  down  the  remaining  amount  at 
the  end. 

N 

4b.  If  the  data  has  been  scaled  up  by  less  than  2 , scale  it  down 

by  1/2  at  each  stage  of  the  IFFT  until  such  time  as  no  more  scaling  is  needed. 

The  remainder  of  the  program  is  reasonably  straightforward.  Gaussian 
and  impulse  noise  are  added  to  a new  s(n)  in  the  interrupt  routine,  as 
soon  as  it  comes  in  from  the  A/D  converter  and  before  it  is  sent  to  Mx . 

Gaussian  noise  is  simulated  by  first  generating  a 9-bit  pseudorandom  number 
with  a flat  distribution.  Eight  bits  are  used  as  a pointer  into  a table  of 
the  midpoints  of  256  equal  areas  over  the  positive  half  of  a Gaussian  dis- 
tribution, and  the  ninth  bit  determines  the  sign.  The  numbers  in  the  table  were 

chosen  such  that  a = .25,  and  values,  when  fetched  from  the  table,  are  multiplied 
by  the  appropriate  constant,  GAUSML,  to  obtain  the  desired  a . Impulse  noise 

is  created  by  adding  to  each  s(n)  a signal  of  height  IMPAMP  for  the  duration  of 
time  IMPDUR,  and  at  a periodic  rate  IMPRTE. 

The  phase  > cj)  , and  the  components  <J>  q,  <f>  and  <fx  , are  all  stored  in  the 
computer  as  fractions,  with  1.0  corresponding  to  360°.  <J>o  is  updated  at  each 
s(n),  by  adding  FRQOFF/IOOOO,  since  an  offset  of  1 Hz  would  be  360°  in  10000 
samples,  to  likewise,  is  updated  by  adding  JITTER/10000  at  each  s(n).  The 
cosine  of  ut  is  multiplied  by  PKT0PK/360  (expressing  degrees  peak  to  peak  as  a 
fraction)  to  obtain  ^ . Phase  hits  are  simulated  by  adding  PHSHIT/360  to  <J> 
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for  those  samples  s (n)  during  the  interval  PHSDUR  and  spaced  by  the  time 


PHSRTE.  The  three  values  cos  ux , cos<j)  and  sin<t>  are  determined  by  table  lookup 
from  Mx  using  the  same  first  quadrant  64  point  cosine  table  as  is  used  for  the 
FFT's  plus  an  additional  interpolation  angle  of  90°/128  whose  cosine  and  sine 
are  stored  in  Md.  From  these  it  is  possible  to  determine  the  cosine  and  sine 
of  any  angle  from  0-360°  to  the  accuracy  of  90/128°  using  the  formulas: 
cos  (a+b)  = cosa  cosb  - sina  sinb 
and  sin  (a+b)  = cosa  sinb  + cosb  sina 

Harmonic  distortion  is  added  just  before  the  final  adjustment  of  the  data 

from  block  floating  to  fixed  point  so  as  to  gain  bits  in  the  cubing  and  squaring 

2 3 

of  the  data.  The  terms  g (n)  and  g (n)  are  added  together,  multiplied  by  the 


parameter  HARM,  and  subtracted  from  g(n)  to  obtain  s(n),  which  when  converted 
back  to  fixed  point,  is  the  final  output  of  the  system. 

A somewhat  intuitive  feel  for  the  effects  of  phase  jitter  and  frequency 
offset  can  be  gained  by  referring  to  Figure  10.  This  figure  shows  a time 
exposure  of  the  impulse  response  of  the  telephone  filter  for  Conus  poor  voice 
with  a)  all  of  the  other  distortions  removed  from  the  system,  b)  only  phase  jitter 
added  and  c)  only  frequency  offset  added.  It  can  be  noted  that  with  a phase  jitter 
interference  the  samples  "jitter"  about  a focal  point,  whereas  with  a frequency 
offset  there  is  a smooth  motion  of  the  wave  as  it  appears  to  continually  pass 
through  and  disappear. 

VII . DVT  Times  and  Space 

Table  3 is  a list  of  the  various  subroutines  in  the  LDVT  in  their  approximate 
order  of  occurrence.  The  table  includes  for  each  subroutine  the  amount  of 
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program  memory  required  and  the  amount  of  time  required  for  its  execution. 

The  total  time  needed  per  frame  is  slightly  more  than  half  of  the  time  available 
and  the  program  uses  up  77%  of  Mp. 

There  are  only  512  locations  in  data  memory,  Md,  of  which  half  is  needed 
for  the  128  point  FFT  real  and  imaginary  buffers.  Fortunately,  no  other  large 
buffers  are  needed  by  the  program,  so  that  an  additional  125  locations  were 
needed  for  the  various  parameters,  temporaries  and  variables,  leaving  one 
fourth  of  Md  unused. 

As  for  outboard  memory,  Mx,  256  locations  are  needed  for  the  filter 
coefficients,  256  for  the  Gaussian  table,  128  for  the  cosine  tables,  and 
4 x 128  for  the  various  speech  buffers,  for  a total  of  1152  or  56%  of 
Mx  (Table  4) . 

VIII.  Summary 

A digital  telephone  channel  simulator  has  been  implemented  on  the  LDVT 
which  operates  well  within  the  real-time  requirements  and  allows  for  comfortable 
margins  in  each  of  the  various  LDVT  memories  (Mp,  Md,  Mx) . Included  in  the 
simulations  are  Gaussian  noise,  impulse  hits,  filter  frequency  response  and 
delay  characteristics,  phase  jitter,  frequency  offset,  phase  hits,  and  nonlinear 
distortion.  The  sources  of  each  of  these  impairments  in  the  telephone  channel 
and  the  techniques  for  measuring  them  have  been  discussed,  as  well  as  the 
limitations  of  the  simulation  compared  to  the  real  situation. 

One  could  easily  set  up  a system  to  send  input  speech  through  a telephone 
channel  simulator  inoneLDVT  into  a vocoder  (LPC,  APC,  etc)  in  another  LDVT 
to  investigate  the  amount  and  type  of  degradations  caused  by  tandeming  a 
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telephone  line  with  a vocoder.  The  telephone  simulation  has  been  implemented 
so  as  to  make  it  convenient  for  a user  to  vary  the  parameters  of  the  system 
online,  making  it  possible  to  determine  which  aspects  of  the  telephone  system 
are  damaging  to  the  vocoder  in  what  way,  and  then  to  develop  the  vocoder  to 
be  more  robust  against  these  impairments. 

Discussions  have  been  evolving  recently  concerning  the  possibility 
of  implementing  a digital  modem  on  the  LDVT.  One  can  imagine  at  some 
future  time  a series  of  LDVT’ s representing  telephone  channels,  modems,  and 
vocoders,  which  can  open  the  door  to  a wide  variety  of  possible  experiments. 
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a)  Modulated  up 
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-w 
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0) 
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b)  Band  pass  filtered  to  exclude  frequencies  below 


c)  Modulated  down  at  receiver 


d)  Low  pass  filtered  to  restore  F(co) 

Fig.  1.  Single  side  band  modulation  demodulation  technique. 
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f(t)  =>  F(w) 




\ 

f(t)  =>  F(u>) 


1 

/ 

/ 

J f(t)  =>  jF(u>) 


/I 


Z(u)  = 1/2  (F (to)  + j F (to) ) 
z(t)  = 1/2  (f(t)  + jf(t)) 

Fig.  3.  Analytic  signal  z(t)  expressed  as  a function  of  real  signal  f(t) 
and  its  Hilbert  transform  f(t).  (See  Text  for  Discussion). 
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Fig.  4.  Single  side  band  signal,  f (t),  expressed  as  function  of  analytic 
signal,  z(t). 
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Fig.  5.  Cosine  wave  at  312  cycles/sec  distorted  by  phase  jitter  at  60  cycles/ 
sec,  30°  peak-to-peak. 
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e 
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Fig.  6.  Nonlinear  voltage  transfer  characteristics  of  two-port. 
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Frequency  (Hz) 

Fig.  7.  C-message  frequency  weighting  function. 
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Fig.  8.  Block  diagram  of  telephone  channel  simulator. 


OUTPUT 

s(t) 
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DATA  FOR  CONUS  POOR  VOICE 


012403 

003720 

PHSRTE 

2000. 

012404 

000001 

PHSDUR 

1 

012405 

000001 

FRQOFF 

1 

012406 

000017 

PKTOPK 

15. 

012407 

000074 

JITTER 

60. 

012410 

000021 

PHSHIT 

17. 

012411 

050475 

HARM 

.6347. 

012412 

001320 

GAUSML 

.022. 

012413 

000000 

IMPRTE 

0 

012414 

000000 

IMP AMP 

.0. 

012415 

000001 

IMPDUR 

1 

012416 

000000 

PKTFRC 

.0. 

012417 

000000 

OFFFRC 

.0. 

PERIOD  BETWEEN  PHASE  HITS  IN  .01.  sec  UNITS 
DURATION  OF  PHASE  HITS  IN  .01.  sec  UNITS 
FREQUENCY  OFFSET  IN  HZ 
DEGREES  PEAK  TO  PEAK  OF  JITTER 
HZ 

DEGREES 

HARMONIC  DISTORTION  FACTOR 
NEW  SIGMA  = GAUSML* . 25 . 

PERIOD  BETWEEN  IMPULSE  HITS  IN  .01.  sec  UNITS 

AMPLITUDE  OF  IMPULSE  HITS 

DURATION  OF  IMPULSE  HITS  IN  .01.  sec  UNITS 

FRACTIONAL  PART  OF  PKTOPK 

FRACTIONAL  PART  OF  FRQOFF 


Fig.  9. 


Parameters  of  synthesizer  reside  in  locations  403-417  of  Md. 
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a)  Unit  Sample  Response 


b)  Same  as  in  (a),  but  with  phase  jitter  at  17°  peak 
to-peak,  60  Hz 


c)  Same  as  in  (a) , but  with  frequency  offset  between 
transmitter  cosine  and  receiver  cosine 
of  1 Hz 


Fig.  10.  Time  exposure  of  unit  sample  response  of  filter  for  conus  poor 
voice  line. 
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TABLE  1. 

Impairments  of  the  Eight  Telephone 
Channels  Simulated 


CHANNEL 

PHASE 

FREQ. 

PHASE 

HARMONIC 

GAUSSIAN 

IMPULSE 

SIMULATED 

HITS 

OFFSET 

JITTER 

DISTORTION 

NOISE 

NOISE 

CONUS 

POOR 

45/15  min 

1 Hz 

15°pk  to  pk 

-20  dbmc 

-40  dbmc 

none 

VOICE 

0 17° 

60  Hz 

CONUS 

MID 

15/15  min 

0 Hz 

11 .6°  pk  to  pk 

-25  dbmc 

-46  dbmc 

none 

VOICE 

017° 

60  Hz 

CONUS 

POOR 

45/15  min 

1 Hz 

14°  pk  to  pk 

-20  dbmc 

-37  dbmc 

none 

DATA 

017° 

60  Hz 

CONUS 

MID 

5/15  min 

0 Hz 

11°  pk  to  pk 

-28  dbmc 

-45  dbmc 

none 

DATA 

017° 

60  Hz 

EUROPEAN 

POOR 

405/15  min 

7 Hz 

35°  pk  to  pk 

-37  dbmc 

-39  dbmc 

225/15  min 

VOICE 

042° 

50  Hz 

74  dbrnc 

EUROPEAN 

MID 

45/15  min 

3.5  Hz 

26°  pk  to  pk 

-43  dbmc 

-45  dbmc 

25/15  min 

VOICE 

032° 

50  Hz 

74  dbrnc 

EUROPEAN 

POOR 

135/15  min 

6 Hz 

35°  pk  to  pk 

-37  dbmc 

-39  dbmc 

225/15  min 

DATA 

042° 

50  Hz 

74  dbrnc 

EUROPEAN 

MID 

135/15  min 

2.7  Hz 

18°  pk  to  pk 

-43  dbmc 

-45  dbmc 

25/15  min 

DATA 

022° 

50  Hz 

74  dbrnc 

TABLE  2. 


Values 

used  for 

the  parameters  of  the 

simulator 

to  correspond  to  the 

data  in 

TABLE  1. 

CONUS 

CONUS 

CONUS 

CONUS 

EUROP. 

EUROP. 

EUROP. 

EUROP 

POOR 

MID 

POOR 

MID 

POOR 

MID 

POOR 

MID 

VOICE 

VOICE 

DATA 

DATA 

VOICE 

VOICE 

DATA 

DATA 

PHSRTE 

(.01  sec  units) 

2000 

6000 

2000 

18000 

222 

2000 

667 

667 

PHSDUR 

(.01  sec  units) 

1* 

1 

1 

1 

1 

1 

1 

1 

FRQOFF 

(Hz) 

1 

0 

1 

0 

7. 

3. 

6. 

2. 

PKTOPK 

(degrees) 

15. 

11. 

14. 

11. 

35. 

26. 

35. 

18. 

JITTER 

(Hz) 

60. 

60. 

60. 

60. 

50. 

50. 

50. 

50. 

PHSHIT 

(degrees) 

17. 

17. 

17. 

17. 

42. 

32. 

42. 

22. 

HARM 

.6348. 

.3455. 

.6348. 

.2563. 

.0928. 

.0464. 

.0928. 

.0464 

GAUSML 

.022. 

.0085. 

.0275. 

.0098. 

.0198. 

.0082. 

.0195. 

.0082 

IMPRTE 

(.01  sec  units) 

0 

0 

0 

0 

400. 

3600. 

400. 

3600. 

IMPAMP 

0 

0 

0 

0 

.11. 

.11. 

.11. 

.11. 

IMPDUR 

(.01  sec  units) 

1* 

1 

1 

1 

1 

1 

1 

1 

PKTFRC 

0 

.6. 

0 

0 

0 

0 

0 

0 

(degrees) 

OFFFRC 

0 

0 

0 

0 

0 

.5. 

0 

.7. 

(Hz) 

PHSDUR  § IMPDUR  are  variables  of  the  system  which,  however,  were  fixed  at  .01.  sec  for  these  simulations  since 
there  was  no  information  given. 


+ 

* 

TABLE  3. 

List  of  Program  Subroutines  and  Time  and  Program  Memory 
Required  in  Telephone  Channel  Simulator 

SUBROUTINE 

FUNCTION 

PROG.  MEM. 

TIME 

GTSPCH 

Fetch  input  speech  from  Mx,  store  in  Md 

29 

.078  msec 

SCLUP 

Scale  data  as  block  left  maximally,  twice 

39 

.422  msec 

FFTFWD 

128  point  forward  FFT 

137 

1 .2  msec 

EOLL 

Bit  reversed  even  odd  separation.  Last  loop  of  256  pt  FFT 

109 

.28  msec 

FLTMUL 

Complex  multiply  by  filter  coefficients 

29 

.28  msec 

SVEBUF 

Save  filtered  spectrum  in  Mx 

19 

.08  msec 

IFFT1 

128  pt.  IFFT  to  obtain  even  numbered  samples 

94 

1.20  msec 

GTBUF 

/v 

Save  f(n),  f(n),  even  numbered,  in  Mx 
Fetch  filtered  spectrum  from  Mx 

44 

.12  msec 

IFFT2 

128  pt.  modified  IFFT  to  obtain  odd  numbered  samples 

-(shared) 

1.20  msec 

FTCEVN 

Fetch  f(n),  f(n),  evens,  back  from  Mx 

14 

.39  msec 

PHASIT 

Compute  phase,  compute  g(n)  = f (n)  cos  <f>  + ?(n)  sin  <f> 
Add  nonlinear  distortion 

121 

.96  msec 

ADJSCL 

Change  s (n)  from  block  floating  to  fixed  point 

29 

.046  msec 

STRSPC 

Store  s(n)  in  Mx 

16 

.039  msec 

ADINT 

Interrupt  routine.  Service  A/D-D/A 
Add  Gaussian  and  Impulse  Noise  to  s(n) 

Update  counters  for  phase  hits  and  impulse  hits 

109 

.46  msec 

TOTAL 


789  = 77% 


6.755  nsec  = 53% 


TABLE  4. 

Buffers  in  Outboard  Memoyr,  Mx,  needed  for  Telephone  Simulation 


BUFFER 

ALLOCATION 

BUFFER  SIZE 

EVIO 

Even  numbered  in-out  samples 

64 

ODIO 

Odd  numbered  in-out  samples 

64 

EVNEW 

Most  recent  64  even  numbered  samples  of  s(n) 

64 

ODNEW 

Most  recent  64  odd  numbered  samples  of  s(n) 

64 

EVOLD 

Previous  64  even  numbered  s (n)  7 

, , . >■  also  as  Mx  Real 

64 

ODOLD 

Previous  64  odd  numbered  s (nj  j 

64 

MXIMAG 

Temp,  storage  of  imaginary  FFT  output 

128 

GAUS 

Midpoints  of  256  equal  areas  of  positive  half  of  Gaussian 

256 

FLTR 

Bit  reversed  filter  coefficients  (128  Real,  128  Imag.) 

256 

SINE 

Sine  of  0 to  it/ 2 for  FFT's  and  for  computation  of  cosine  4> 

64 

RVSINE 

Bit  reversed  sine  of  0 to  tt/2  for  last  stage  of  forward  FFT 

64 

TOTAL 

1152  = 56% 
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